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Abstract 

We explore the non-equilibrium dynamics of two coupled zig-zag chains 
of trapped ions in a double well potential. Following a quench of the 
potential barrier between both wells, the induced coupling between both 
chains due to the long-range interaction of the ions leads to their complete 
melting. The resulting dynamics is however not exclusively irregular but 
leads to phases of motion during which various ordered structures appear 
with ions arranged in arcs, lines and crosses. We quantify the emerging 
order by introducing a suitable measure and complement our analysis 
of the ion dynamics using a normal mode analysis showing a decisive 
population transfer between only a few distinguished modes. 


1 Introduction 

Cooled ions in traps form a clean and highly versatile setup for exploring struc¬ 
ture formation with long-range interacting particles, both in equilibrium and 
non-equilibrium. In equilibrium, possible structures include small and large 
ion crystals Elm [7] possessing various internal ordering such as concentric 
rings (2D), shells (3D) [3J[ni[l] and string-of-disks configurations [52], and even 
two-component Coulomb bicrystals m- At the crossover from one to higher 
dimensions trapped ions can also form zig-zag configurations, a structure that 
attracts particular attention in the recent literature [11130112]. Specifically, vary¬ 
ing the geometry of an anisotropic harmonic trap allows for a second-order phase 
transition from a linear to a zig-zag structure, which can be either ideal or in¬ 
volve topological defects [iniiMiiis]- Out of equilibrium, the recent literature 
predicts an equally rich variety of possible ionic structures including spatiotem- 
poral patterns in laser-driven microtraps [24] and periodic lattices [5ni[57], but 
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also interaction induced current reversals of the transport direction based 
on structure formation in the phase space. Much of this research on structure 
formation with trapped ions roots in the admirable advancements of the con¬ 
trollability of ions in recent years. This ranges from the quickly progressing 
miniaturization of ion traps and lab on chip technologies nasD] via the advent 
of optical trapping techniques [3S] to the discovery of multi-segmented Paul and 
Penning traps [37113S]. The latter example in particular allows for more and 
more complex but still controllable arrangements of long-range interacting par¬ 
ticles as required e.g. for quantum information processing UmiSlISeKIS]. The 
above advancements allow and evoke a new type of question: How do individual 
ionic structures respond if we couple them to each other? Consider for example 
a segmented ion trap with two wells, both loaded with an individual zig-zag con¬ 
figuration, separated from each other by a potential barrier between the wells. 
Let us now quench the barrier to a lower value, which increases the coupling 
between the individual chains: Are the only two (expected) alternatives for the 
dynamics of the zig-zag chains that they either deform only slightly and re¬ 
spond with small oscillations to the increased coupling or that we observe their 
complete melting resulting in irregular oscillations of all ions? This is precisely 
the problem we want to investigate in the present work. To explore the above 
problem we develop a minimal model based on a two dimensional double well 
potential that allows for zig-zag configurations in both wells whose geometries 
resemble the well-known zig-zag states in anisotropic single well traps. Remark¬ 
ably, following the quench of the barrier we observe that the melting process 
does not simply lead to irregular oscillations but to a complex non-equilibrium 
dynamics constituted of different phases of motion. Phases of irregular oscilla¬ 
tory motion are interrupted by motional phases which exhibit transient ordered 
configurations. Although nonlinear dynamics governs the motion of the cou¬ 
pled ion chains, we employ a normal mode analysis showing that the population 
of the corresponding linear eigenvectors is not arbitrarily distributed over the 
whole band of modes as one would expect for e.g. a chaotic system. Instead, 
during most phases of the time evolution only a few eigenvectors are strongly 
populated and nonlinear effects show up in form of a decisive and quite sud¬ 
den transfer of energy among the different modes. Our work is organized as 
follows. Section II explains our setup and the preparation of the ground state 
configuration. Section III provides our main results followed by a normal mode 
and population analysis of the dynamics. We summarize our findings and their 
interpretation in section IV. 


2 Setup, Hamiltonian and ground state config¬ 
uration 

We consider N ions, described as classical point particles with mass m and 
charge Q, confined in radial direction {x,y) to a linear quadrupole Paul trap 
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Figure 1: (a) Cartoon of the double zig-zag equilibrium configuration in the double 
well potential used as the initial configuration. (b,c) Magnifications of the part of the 
ionic configuration which links between the two chains, i.e. in the barrier region. It 
highlights the difference between the ground state configuration (b) and the energeti¬ 
cally next higher equilibrium configuration (c) which we call the mirror configuration, 
(d) The red (solid) line shows the double well potential before the quench and the blue 
(dashed) line after the quench. A lowered barrier enhances the coupling between the 
two ion chains (arbitrary units throughout). 




and to a double well potential (segmented trap) in z—direction. 

= ^{cx^ + cy‘^) + ^cos(wrft)(ca;^ - cy^). (1) 

with C/dc and C/j-f being the applied constant and the rf-voltage; Wj-f is the (ra¬ 
dio) frequency, c is a geometrical parameter of the trap. The geometrical param¬ 
eter is for both directions x and y equal in radially symmetric traps but in planar 
traps [m [32] the geometrical parameter for both directions can strongly differ. 
The ion dynamics in the radio-frequency trap is composed of the so-called micro 
motion, and a comparatively slow averaged motion taking place in an effective 
harmonic potential [T2| V{x,y) = ^(w^x^ -ha;^?/^). Here, Ux = — q‘^/2 

and ujy = -I- q‘^/2 are the effective trapping frequencies with a = ^^^ 2 “ c 

and q = being dimensionless parameters. For the confinement in z- 

direction we assume the following phenomenological double well potential [3H] . 
with wells centered at « ±zo and separated from each other by a barrier of 

^Note that formally, the Laplace equation does not allow to combine an ideal linear Paul 
trap with a double well potential in axial direction. However, experimental traps are not ideal 
and reference m could indeed realize the combination of an approximately linear Paul trap 
and a Mexican hat like potential which justifies to consider a combination of a double well 
potential with a linear Paul trap. 
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height ^ IjC (see Fig.[^. 
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This potential quantitatively resembles the shape of individual harmonic wells 
around ±zo up to terms proportional to . 

Specifically, for a given sufficiently high barrier this allows us to prepare zig¬ 
zag chains in each of the two wells which are the energetically lowest equilibrium 
configuration of the double well which we call the ground state configuration 
in the following Note that finding the many-particle minimum of a many ion 
system is generally a highly nontrivial task; hence the present choice of the 
double well potential is a crucial step to allow for a numerical study of the 
dynamics of coupled ion chains. After preparing this configuration, i.e. its 
numerical determination, our strategy will be to ramp down the barrier height by 
a certain amount which corresponds to a quench of the quantity C. Subsequently 
the resulting dynamics of the now strongly coupled ion chains will be explored. 
To understand the complex dynamics of coupled many-ion structures in non¬ 
equilibrium it is crucial to simplify our model. First, since we are interested 
in the dynamics on large scales we neglect the micromotion. Second, we focus 
on a two-dimensional description, which simplifies the visualization of the ionic 
structures and their analysis but does not change qualitatively the resulting 
dynamics and phenomenology of the structure forming processes. Specifically, 
we choose parameters a = uJx/uJz ^ S.3 where the ground state configuration is a 
planar (2D) zig-zag structure in the a;—z-plane far from the transition to a helical 
(3D) zig-zag chain. The transition from the 2D structure to the 3D helical zig¬ 
zag chain occurs at a « 4.5 , whereas the transition to a linear ID line-structure 
occurs at a « 10 (for ujx = ujy)- Thus, although ramping the barrier will 
generally produce both in-plane and out-of-plane fluctuations, the latter ones 
are generically small. A strong confinement in y-direction (wy/w^ > 10) prevents 
their amplification in the course of the consecutive dynamics; hence the ions stay 
close to the x — z plane. Indeed the impact of small out-of plane fluctuations on 
the intra-plane dynamics is quadratically suppressed with the distance of the 
ions perpendicular to the x — z plane and consequently a projection of the ion 
dynamics in 3D looks very similar to a direct 2D description. Hence, we focus on 
a two-dimensional minimal model given by the two-dimensional Hamiltonian: 


n 
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(3) 


with V{x) = miJ^x^ jl, = {xi,Zi), 

^The only force which couples the dynamics of ions perpendicular to the x — z plane 
to the intra-plane dynamics is the Coulomb coupling. Consider two ions in distance L = 
.^L| -|- where Ly is the projection of their distance onto the x — z plane and is their 

distance perpendicular to this plane. Then Fy = T’cos[tan(Lj^/L||)] ps T’[1 — (Lx/L||)^/ 2] for 
L|| ^ showing that forces produced by small out-of plane oscillations onto the intra-plane 
dynamics are quadratically suppressed. 
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Tij = ^ {xi — Xj)"^ + (zi — ZjY and p = {px,Pz). Introducing a rescaled time 
tu = ^luJz and space units Xu = K = [(^^/(dTreoma;^)]^^^ and defining 


t = uJzt; X = —; 


K' 


K’ 


K’ 


C* = 


c 






a = 


(4) 


we integrate the resulting Newtonian equations of motion. For the sake of 
convenience we will drop the stars in the following. Finding the ground state 
configuration is a 2N dimensional optimization problem which is, for large iV, 
numerically highly demanding, mainly because root finding algorithms converge, 
for most initial configurations, to some excited equilibrium configuration. To 
avoid this, we need a good initial guess for the positions of the ions which 
converges to the ground state configuration. Here, the quantitative similarity 
of both wells of the chosen double well with individual harmonic traps comes 
into play: it allows us to exploit the existing knowledge on a single zig-zag ionic 
chain in the literature and hence to numerically efficiently determine the double 
zig-zag ground state configuration within standard multidimensional rootfinding 
which is visualized in Fig. Note that the present setup also allows for 
a ’first excited’ equilibrium configuration which is very similar to the ground 
state configuration (Fig. Sb, c)) and whose energy exceeds that of the ground 
state configuration by a factor of less than 10“^°. We call this first excited 
equilibrium configuration the mirror configuration. 


3 Results and Analysis 

We now use the planar double-zig-zag ground state configuration as our initial 
configuration and assume that the initial velocities of all ions vanish which 
can be achieved approximately via modern cooling techniques in experiments 
with trapped ions [HI m 0 [HI mi- We excite this initial configuration by 
performing a sudden quench of the barrier height between the two wells. This 
allows in the course of the dynamics the zig-zag structures in both wells to move 
closer together, thereby effectively increasing their coupling. Here we quench, 
exemplary, from C = 1.26-10“^ to C = 6.3. This quench leads to an energy 
excess of less than two percent with respect to the ground state configuration in 
the C = 6.3 potential which is still a double zig-zag configuration. To explore 
the time evolution of the double zig-zag configuration after the quench we use 
an implicit Gaussian 4th order Runge-Kutta algorithm m for integrating the 
Newtonian equations of motion. For very short times we mainly observe a center 
of mass motion of both zig-zag structures towards each other, accompanied by 
small oscillations of the individual ions. At t ~ 1.5 the two chains begin to 
melt in the center of the system which continues until the zig-zag structures 
are completely destroyed. Strikingly, as a result of the melting process, we 
do not only obtain phases of irregular oscillations of the ions in the trap (see 
for example Fig. [^at t = 18.91), but also many phases of transient structures 
featuring an unexpected degree of order: we observe the formation of ions in 
lines, arcs and cross like structures and partial revivals of the dynamics (see 
Fig.§. As a result, coupled zig-zag chains not only exist as a stable ground 
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and excited equilibrium configuration but, once exposed to perturbations such as 
our quench, experience a complex reordering process with both irregular phases 
of motion and ordered configurations as quantified in Fig. [^a). It turns out 
that the shape of the specific configurations we observe is sensitive to the initial 
configuration: Choosing the mirror configuration instead of the ground state 
configuration as our new initial configuration and performing the same quench 
of the barrier as before leads to a second set of ordered structures, including arc 
like configurations and elliptical arrangements of ions (see Fig.j^b). 
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Figure 2: Snapshots of the time evolution of = 44 ions in the double well for 
an initial configuration given by the ground state configuration (left column) 
and the mirror configuration (right column) before the quench. Parameters (in 
scaled units): a = 8.25; Zq = 5/2^/^; C is changed from 0.0126 to 6.26. 


In the following we introduce a measure based on the idea of Voronoi dia¬ 
grams m to quantify the degree of order underlying the structures shown in 
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Fig-i Specifically, when drawing a circle around each ion, whose diameter is 
given by the distance to its nearest neighbor, our measure is defined as the sum 
of the areas of the respective circles around all = 44 ions, i.e. 

Mt) = ^ + « W - (5) 

i 

To evaluate this measure we rescaled the z-coordinates of all ions z* = z ■ k 
with k = 0.04, such that the average x-distance of adjacent ions equals their 
average distance in z-direction. The result of A{t) is shown in Fig. a) and 
features pronounced oscillations: all structures which appear to be ordered from 
a visual impression and show ions arranged in arcs, lines or crosses (Fig. lead 
to distinct minima in A, whereas irregular structures (Fig. [^d (left and right), 
f (left) and g (right)) generate larger values of A. Clearly A is a measure for 
the clustering of ions on one-dimensional manifolds such as straight lines, arcs 
or crosses leading to low values of A, whereas acquires larger values if the ions 
tend to cover a two dimensional areaj^ We now analyze the observed emergence 



Figure 3: Time evolution of the measure A as defined in Eq. (5). Red squares 
show times at which snapshots in Fig. 2 were taken for the ground configura- 
tion(left column). Insets illustrate the definition of A as the sum of the areas of 
all shown circles. Upper and lower insets illustrate how A accounts for typical 
irregular and regular configurations respectively. 

of transiently ordered ion structures originating from the melted double zig-zag 
structures. Clearly, if the initial quench of the barrier had induced only small 

®Note, that A would be also small if all ions would cluster, which is however excluded from 
the present situation due to the chain character of their arrangement. 
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oscillations around the initial configuration of the ions the dynamics could be 
described to a good approximation by the linearized equations of motion. In 
terms of normal modes |13j , it would then consist of a superposition of harmonic 
oscillations with preserved oscillation amplitudes. 

Here, the ion dynamics does not consist of small amplitude oscillations (the 
initial configuration dissolves completely) but is essentially nonlinear. Gener¬ 
ally, nonlinearities couple different normal modes and transfer energy between 
them resulting in varying normal mode amplitudes. Here, it turns out how¬ 
ever, that a normal mode analysis is still useful: We find that during certain 
phases of the ion dynamics only a weak energy transfer between different normal 
modes takes place, indicating that nonlinearities have only a weak impact on 
the ion dynamics during these phases. These phases are interrupted by short 
transition ’events’, where significant population is transferred between different 
eigenvectors. Specifically, we use standard techniques to numerically calculate 
the eigenfrequencies {uJi} and the associated eigenvectors {E^} around the new 
many-ion equilibrium configuration of the double well following the quench. As 
the eigenvectors form a basis of the Wdimensional configuration space (here 
sorted by their eigenfrequencies) we can now represent position and velocity of 
each ion at a fixed time t in this basis: 


(ri,..,rjv) (t) = (r®"‘,..,r^)-P^p,(t)Ei; 



( 6 ) 


Here, is the position of the m-th ion in the ground state configuration (mirror 
configuration) of the double well following the quench and pf{t) {k‘f{t)) describe 
the population of the linear eigenvectors E^. In the linear regime the population 
coefficients reduce to pi{t) = Pi cos{uJi * t) and ki[t) = Ki sin(a;i * t) with Pi, Ki 
being constant amplitudes. 

We now analyze the population dynamics of the eigenvectors E^ as induced 
by nonlinear effects. The initial quench mainly populates the vectors E2, E4, Eg 
(Fig. 1^ at t = 0). These belong to the slowest modes whose population is 
allowed by symmetry (they possess the same symmetry as the initial configu¬ 
ration). Fig. ([^) shows that a first significant energy transfer to higher modes 
takes place at t ~ 1.5, which is the time where the zig-zag configurations start 
to melt, and then spreads to more and more other modes. Notably, the energy 
does not spread in an arbitrary manner over a whole band of modes which we 
would expect for a purely irregular or strongly chaotic system, but many of the 
eigenvectors either keep their population for relatively long times (Fig. |^) or 
exchange them only pairwise (Fig. [^. This indicates that the many ion dynam¬ 
ics is not irregular but possesses a high degree of order. We now study in more 
detail how the population transfer between individual linear normal modes takes 
place. Therefore, we populate only one (or several) of the linear eigenvectors 
as our initial configuration and explore the resulting dynamics. Populating for 
example only vector E 2 which has the largest population after the quench, i.e. 
choosing (ri,.., rjv) (t = 0 ) = (r®‘^,.., r^) + = 0)E2 leads to a permanent 

’breathing’ oscillation of all ions in the left well and an opposite motion of the 
ions in the right well. Accordingly, we observe almost no transfer of popula¬ 
tion to other eigenvectors and the considered single mode excitation leads to 



Figure 4: Time-dependent population of linear eigenvectors defined as Pi{t) = 
{di{t) cos(uJit))^ (color) for different initial populations: a) created by the quench 
for the ground state configuration as initial configuration, b) as in a) but af¬ 
ter setting all but the three largest initial populations {P 2 ,P 4 ,P 6 ) to zero, c) 
Population of only one mode with P 2 &{t = 0) = 0.11524. Parameters like in 
Fig.H 


persistent small and periodic oscillations rather than inducing a melting of the 
double zig-zag. If we populate, e.g., the three eigenvectors which attract the 
strongest population by the quench of the potential barrier, we obtain a similar 
result: Only very few energy is transferred to other modes (Figj^). There are 
however eigenvectors belonging to higher modes, whose initial population, even 
if it is small, induce a hierarchical spreading of population between different 
eigenvectors. In particular, if we weakly populate vector E 26 above P 2 Q{i = 0) 


9 









-0.4 


-0.5- 

0 5 10 15 20 25 30 


t 


Figure 5: Magnification of the population dynamics of (i 4 o(t) cos(a; 4 ot) and 
d 4 i(t) cos(w 4 it) as defined in eqj^ 


we obtain a population dynamics which reproduces large parts of the occupation 
dynamics as induced by the quench of the potential in the full problem (Fig.|^). 
Already a weak population of this vector does therefore lead to significant non¬ 
linear effects strongly mixing the linear eigenmodes. As a result, it depends 
on the details of how we excite the double zig-zag ground state configuration 
whether it responds with small oscillations or with a complete melting and the 
dynamical reconfiguration of the ions into ordered (transient) structures. Here, 
we neglected micromotion for simplicity. Accounting for it would however influ¬ 
ence the occupation of eigenvectors |20j , but presumably only on relatively long 
time scales. Accordingly, it could be interesting to explore if micromotion leads 
to a diffusion-like spreading of population over many eigenvectors and therefore 
generates a slow crossover from the partly regular dynamics quantihed above to 
an increasingly irregular ion dynamics on long timescales. 

Let us finally address the experimental realization of our setup employing state 
of the art ion technology. Typical experimental parameters for segmented Paul 
traps are Wf.//27r = 4.2 — 50MHz, Urf = 8 — 350V with applied dc voltages in 
axial direction of up to lOV [mill]. Depending on the ion species and trap de¬ 
sign this results in radial direction in a;/27r = 1 — 5MHz and in axial direction in 
= 0 — 5MHz. For the dynamics only the frequency ratio a = — = 8.25 is 
of importance (e.g. a;/27r = 4.5MHz and = 0.545MHz). The parameters 

zg and C depend on the dc voltage and the trap geometry and are in the range 
of 30 pm for the well position Zg and up to 300 pm^ • MHz^ for C. The ion 
configurations in the course of the dynamics could be detected by fluorescence 
light with a CCD camera [211 [H]. 
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4 Conclusions 


While structure formation with trapped ions is already a versatile and active 
research area, recent experimental progress suggests the coupling of individual 
ionic structures in multi-segmented traps. The present work elaborates, on a 
minimal model that allows for coupled zig-zag ion chains that respond in a highly 
selective way when we excite them: Depending on the details of the excitation 
we either obtain regular oscillations around the equilibrium configuration or a 
complete melting followed by a sequence of structures with ions arranged in 
lines, arcs and cross like formations. Our model can be extended to explore the 
dynamics of whole arrays of coupled ionic structures in multisegmented traps. 
The selective response of coupled zig-zag chains to weak excitations promises an 
energy transport (or sound propagation) along the trap which strongly depends 
on the initial signal. Particularly in cases where the ions respond with small 
oscillations we expect a much weaker energy transport than in cases of complete 
melting where nonlinear effects couple linear eigenmodes and convert Coulomb 
energy to intersite from one site of the double well to the other one motion of 
ions. A further application of coupled ionic structures in multiwell traps could 
exploit the current interest in topological defects (kinks) in zig-zag structures 
to study their collision and transferability between different individual zig-zag 
chains. 
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